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Axial oscillations relevant to the r-mode instability are studied with slow rotation formalism in 
general relativity. The approximate equation governing the oscillations is derived with second-order 
rotational corrections. The equation contains an effective 'viscosity-like' term, which originates from 
coupling to the polar g-mode displacements. The term plays a crucial role on the resonance point, 
where the disturbance on the rotating stars satisfies a certain condition at the lowest order equation. 
The effect is significant for newly born hot neutron stars, which are expected to be subject to the 
gravitational radiation driven instability of the r-mode. 



^ - I. INTRODUCTION 

o ; 

\ The surprising discovery of the r-mode instability in rotating stars has inspired the study of the axial oscilla- 

I tions |p]-|l9|. The physical mechanism of the instability is the same as that in the polar modes, so-called the radiation 

?H reaction instability found by Chandrasekhar |^ , and Friedman and Schutz . The r-mode oscillations seem to be 
more important since they are unstable on an inviscid rotating fluid even for small angular velocity. This instability 

' can explain the spin-down process of newly born neutron stars which rotate with nearly Kepler frequency. The grav- 

QQ , itational waves associated with the unstable r-mode oscillations may be promising detectable sources on the ground 

' based laser interferometric detectors. It was proposed that the unstable mode might also play a key role on the 

\ spin of the accreting white dwarfs |l^Jl^. However, Lindblom ||l5[ showed that the possibility is unlikely realized. 

T-H It was also proposed that the r-mode instability could provide the loss of angular momentum from the accretion 

^ ' disk by the gravitational waves to halt the spin-up in the low mass X-ray binaries |13 1^. The proposed conclusion 

o 



crucially depends on the poorly understood dissipation mechanism |17| . In this way, the r-mode oscillations enrich 
the astrophysical implications. Recent review of this subject is given by Friedman and Lockitch [ p^ . 

Most of the studies are however limited to idealized situation. Some effects are added to the simplified models to 
' examine the validity. For example, RezzoUa, Lamb and Shapiro JTsf suggested that the magnetic field of a neutron 
' star is wound up during the non-linear growth of the unstable mode, and that the energy is not transfered to the 
' gravitational radiation so much. However, their calculation is not self-consistent non-linear one, so that the magnetic 
\ effect on the r-mode is not conclusive at moment. The relativistic effects are also important for the oscillations in 
Q^,. neutron stars. The relativistic factor is of order 0.2, so that the frequency could slightly shift from the Newtonian 
^ calculation. More important effects of the general relativity are gravitational wave and frame dragging. Each of them 
(3jr), leads to qualitatively a different result. Some authors |^,^ already calculated the r- modes in general relativity. The 
L| ■ frame dragging effect is regarded as a kind of differential rotation The rotational effects are however limited to 
\ the lowest order. Mathematically, the treatment is insufficient, since the modes are degenerate at the order. It is 
• necessary to include the higher order rotational corrections. Such a task will be significantly complicated when the 
5h ' rotation and relativistic effects are simultaneously considered. 

. . One simplification for the non-rotating case is realized as the decomposition of the spherical harmonic function. 

Each oscillation mode can be specified for each index l^m. Furthermore, the polar and axial modes are definitely 
determined by an appropriate combination of the harmonic functions. However, the functions should be mixed in 
the presence of the rotation. Considering the slow rotation, the coupling is weak, so that the entangled range can be 
restricted. It is not known a priori to specify the axial oscillations on the rotating stars by a few spherical harmonic 
functions. Indeed, Lockitch and Friedman [lO| calculated the normal mode by the sum of infinite number of the 
spherical harmonics indices. We will consider a different approach in this paper. Our treatment is suitable for the 
initial-value problem. We do not consider a single Fourier mode e~*'^* with respect to time. Suppose that the initial 
disturbance at < = is produced with a certain symmetry, which can be assumed to be specified by a few number of 
spherical harmonic functions. What happens in the subsequent evolution? If the oscillation possesses the symmetry, 
the oscillation preserves the symmetry. Otherwise, new patterns with different spherical harmonics indices are induced. 
It is clear that the truncated approximation to the finite number of the spherical harmonics becomes worse for large t 
in general. Therefore, our method is constructing the pulsation equation adequate for small t. It is difficult to address 
the valid domain of t beforehand. For example, the r-mode oscillation can be well described by a few number of the 
spherical harmonics indices for the oscillation in the uniform density with the Cowling approximation [0j2^ . 

Present authors applied the method to the axial oscillation in a rotating relativistic star with the Cowling approx- 
imation 0, where the metric perturbations were neglected. They took account of the rotational correction up to 
third order to examine the oscillation equation. They pointed out the importance of the frame-dragging effect, which 
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causes different property unlike the Newtonian case. Thus, the perturbative approach proves useful in providing a 
physical understanding of many processes. In this paper, we extend the approximation scheme to include the metric 
perturbations. The remainder of the paper is organized as follows. In Sec. II, we formulate the perturbation scheme to 
solve the linearized Einstein equations. We employ the slow rotation approximation, i.e., angular velocity is assumed 
to be small expansion parameter. We look for the solution in which axial-led functions are dominant, where the 
axial-led functions mean the functions describing axial modes in the absence of the rotation. We use the terminology 
polar-led in the same way. In our scheme, the lowest order equations are determined only by the axial-led functions. 
They were already derived elsewhere |^,||,^ 23 1 but are reviewed in Sec. III. In Sec. IV, first-order corrections to the 
polar-led functions are shown. In Sec. V, second-order correction terms are added to the same equations in the form 
as the lowest ones. In Sec. VI, concluding remarks are given. Throughout this paper, we work in the geometrical 
units of c = G = 1. 



II. PERTURBATION SCHEME 



We consider a rotating star with uniform angular velocity f2 0(e), where e is a small rotational parameter. The 
metric and fluid quantities describing the equilibrium state can be calculated by the slow rotation formalism |^,^ . 
They are summarized in the Appendix A. We next investigate the perturbations from the state. The metric pertur- 
bations can be described by six functions. Working in the Regge- Wheeler gauge [|6|, the perturbations are expressed 

as 
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where Yim = ^/m(^, 0) represents spherical harmonics. The symbol "sym" indicates that the missing components of 
hy^i, are to be found from the symmetry h^J^l, ~ /i^^. The angular part is expanded with an appropriate combination 
of the harmonics. In a same way, the fluid perturbations are described by five functions. They are the pressure 
perturbation density perturbation 8p and three components of the 4-velocity (Sur^ Sug, Su^). The component 6ut 
can be determined by the condition, 5u^u^ — 0. These perturbed quantities are also expanded as 
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These eleven functions are determined by ten components of the linearized Einstein field equations, 
and the adiabatic condition for the perturbations. 



Sp + ^-Vp = 
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{6p + C-Vp), 
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where F is the adiabatic index and ^ is the Lagrange displacement. 

Now we will solve the pulsation equations by the expansion of the spherical harmonics. In the spherically symmetric 
star, the equations are decoupled for each harmonic index (I, m). The perturbations can also be decoupled into the axial 
and polar perturbations. They are respectively described by the axial functions Aim = {Uim-, ho^irm hi^im)-, and the 
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polar functions Vim = {5pim,Spim, Rim, Vim, Ho^ijn,Hi,irn, H2,im,Kijn)- In the presence of rotation, the perturbations 
are described by the mixed state of them. From now on, we call these functions as the axial-led ones for Aim and 
the polar- led ones for Vim- Since the slow rotation is associated with the perturbation with ^ = 1, the formal relation 
between Aim and Vim in the Eqs.(0)-(^) can schematically be expressed as 

= [Aim] +£X [Vl±lm] X [Alm,Al (9) 

= [Vim] +SX [Al±lm] X [Vlm,Vl±2m] + ' ' ' , (10) 

where the symbol £ denotes some functions of order e, and the square bracket formally represents the relation among 
the perturbation functions therein. These selection rules follow from the addition of angular momenta. We moreover 
assume that the axial-led functions are dominant in the slowly rotating star, i.e.. Aim S> Vim- This assumption is not 
valid for the some cases. If the star and its perturbations obey the same one-parameter equation of state, r-modes 
and g- modes are coupled at the lowest order in general. They are respectively described by Vim and Aim, which 
have zero frequency in the spherical isentropic star. They form hybrid modes in the rotating star ||l0| . However, the 
r-modes are discriminated from the g-modes for non-isentropic stars, since the g-mode frequencies are non-zero in the 
spherical limit. In that case, we may use the assumption Aim ^ Vim to solve the r-mode oscillations, and expand as 

Aim - A'^l + e'A^l + ■--, Vim - eivi^^ + e^v['2 + •••)• (11) 

Substituting these functions into Eqs.(p[)-(|l0|), and comparing each order of e, we have the following equations of 
e"(n = 0,l,2), 

(12) 
(13) 

-£'x[AZa''^22J (14) 

(15) 

We have here assumed that the perturbation in the lowest order is described by a single component of spherical 
harmonic, that is, A^j^l^ = 0, for V ^ I. The polar-led functions in Eq.(^4|) are eliminated by Eq.(13). Equation 
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( p^ ) represents the axial oscillation in the lowest order. Equation ( [Tsl ) is the second-order form of it, and the term 

S'^ X a'iH can be regarded as the rotational corrections. The method to solve the equations is straightforward. The 
first-order equations are solved by the axial-led functions. The polar-led functions are expressed using them. We 
have the second-order equations with the corrections expressed by the axial- led functions in the lowest order. These 
equations are successively solved in the following sections. In the actual calculations, we also assume that the time 
variation of the oscillation is slow and proportional to Q, i.e., 9t ~ f2 ~ 0(e). This is true in the r-mode frequency, 
dt ~ il-2/l{l + l))mn. 

III. LOWEST-ORDER CALCULATION 

In this section, we review the equations governing the axial oscillations at the lowest order. The radial functions are 
decoupled for each spherical harmonic index. We denote it as L = {l,m) for the abbreviation. The relevant functions 
are /io,Li^i,l and Ul- They are calculated from three components, essentially {t(j)) (r(j)) and (9(j)) components, of the 

Einstein equations. We define a function <i>^"'^'' as 

= ^, (16) 



where the superscript denotes the lowest order term in Eq. (pi]) . The relation between the metric functions is given 
by 
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where dr = dt -\- imfl denotes a time-derivative in a co-rotating frame and a prime denotes a derivation with respect 
to r. The axial velocity function is expressed by two ways: 
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Eliminating C/|^^ in Eqs.@-(|l9l), we have a second-order differential equation for <i>]^ , 
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The linear operator C is defined above. This equation is reduced to a singular eigen-value problem 
the time dependence for the mode as the form e"'*^*. The equation can be written as 



1 
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where 
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1(1 + I) 
2m 



(tr — mfi), 
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by assuming 

(24) 
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There is a singular point 7'o in Eq.(|2J) unless (/(rp) = 0, corresponding to the real value of = n7(ro). It is evident that 
the singularity originates from the mismatch in Eq. (|2^) , i.e., the first term vanishes whereas the second never does. 
When the first term, which is formally of first order, is small enough, then higher order corrections become important. 
This situation is very similar to the inviscid shear flows. When the viscosity is small enough, the stability is almost 
determined by the Rayleigh equation, the perturbation equation for the inviscid theory. The Rayleigh equation has 
critical points for some mean fluid. The viscous corrections should be included to determine the behavior near the 
neighborhood of the critical points. Therefore, the equation should be replaced by the Orr-Sommerfeld equation 
derived from the Navier-Stokes equation. 

As we will show in the subsequent sections, the function <i>^^'' can also affect the polar-led functions. When one 

considers the equations of the next order, additional terms appear in the form (|2^). The terms depend on 

different aspect of the background flow, and play an important role as the viscosity terms. 



IV. FIRST-ORDER CORRECTIONS IN POLAR-LED FUNCTIONS 

In this section, we will derive the equations governing the polar-led functions, X^^), (5p(i), 

_R(^) and V'^'^\ As shown in Sec. II, the functions with {I ± 1, m) are coupled with the axial-led functions with (Z, m). 
We will shorten the subscript of the spherical harmonic index as e.g., ip^^ = ^vi^i^-^- These eight functions are 
determined by seven components of the linearized Einstein field equations and one thermodynamical relation. The 
calculations are straightforward, but results are sometimes messy. The pressure ^p^"* and density perturbations 
are expressed by 

47r5pW ^ 2^(po + Po)(if^'i + 2T±l]r2e--$W) ^ 2S±wUi'\ (27) 
4^5p« = _47r^(ij(;i + 2T±r!r2e-''$W) 

-4^±^ (e^/^n7[/«)' + 2r±-^(.nr2e-)'C/«, (28) 
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where 



S+ = j^Q+, 5_ = i±lQ_, T+^IQ+, T_^-{l + l)Q_, (29) 
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These quantities are expressed by the axial- led functions (<&^"^\ u'j^^) and H^^^. There is another relation among Hq^, 

K^:'' and in the field equations. Eliminating Sp^^'^ by Eq.(^8|), we have a second-order differential equation for 
the metric perturbations. The equation can be regarded as relativistic version of the Poisson equation, V'^Scf) = AnSp, 
for the gravitational potential Scj) and the density perturbation Sp. In the relativistic case, the Newtonian potential is 
replaced by -ffg^'' or K^^K The second-order differential equation is explicitly given by 
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When the axial-led function at the lowest order <i>^"'^^ is given, the functions and k'^^ are solved with appropriate 

boundary conditions. In a similar way, we can solve the other polar-led functions, i?^^]-, tl2^±, -R±^ and V^^^ by 

{^'^j} .u'^i}) and (iJg^j., if^^). The expressions for these four polar- led functions are omitted here, since they are 
eliminated in the following calculations, and never appear in the final results. However, here is a comment on using 
the adiabatic condition Eq.(p[). The time derivative of it can be written as 
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This thermodynamical relation determines the function i?^'' unless the Schwarzschild discriminant A vanishes. Oth- 
erwise, we have one constraint for the function u'j^^ through Eqs.(27)-(p8|), and the function should be specified 
in another way. (For example, see the method by [Q in the Cowling approximation.) The mathematical drawback 
for the isentropic case ^ = is related with the coupling of the g- modes. Both r- modes and g- modes are degenerate 
to zero frequency in the non-rotating star, and hence a particular treatment is necessary [ p^ . From now on, we will 
consider the case only. 



V. INCLUDING SECOND-ORDER CORRECTIONS 



So far we have considered ten components of the Einstein equations and one thermodynamical relation. They were 
limited to the lowest-order form with respect to the rotational parameter. In this section, we consider how the next 
order terms modify the equation (23) derived in Sec. III. The relevant equations for this purpose are three components 
of the Einstein equations, i.e., {t(f>) {rip) and (00) components, which are used in the leading order equation in Sec. III. 

These equations contain the relations among the axial-led functions of second-order ^, h\ 'j^, Uj^ and the polar- led 
functions calculated in the previous section. We follow the same procedure as done in Sec. III. Defining the function 



^^^^ = h'^^'j^/r'^ and eliminating hYi^ and 11}^', we eventually have the equation governing the axial oscillations. It can 
be written in the following form 
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where the operator £ is defined in Eq.(|2 
They consist of several terms as 



and the right hand side means the second-order rotational corrections. 
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The explicit forms of the coefhcients a^. Pi, Ai and Bi are given in Appendix B. They are expressed by the quantities 
determined by the stellar model in the equilibrium. Since the term contains higher order derivatives, we explicitly 
show below: 
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a function $i = + e^^^^^ satisfies the following equation, which is correct up to 0(£^), 
C[<Pl] = V[<i>L, /iix, Ul] + g[Ho^±,K±]. (44) 



The quantities without the superscript satisfy the same relations as in the leading order, i.e., Eqs.(pT[), (p^), ( |3l| ) 
and {p2[), which are adequate approximation to this order. This equation is of course reduced to the leading order 
equation (p3|), when the second-order rotational effects and the coupling to the polar modes are neglected. 

The first term in Vq contains fourth derivative of with respect to r, since Ul can be expressed by the second 
derivative of $l as shown in Eq.([l9|). The highest derivative term of $l with respect to r is therefore given by 
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Neglecting Q and V except 2?* in Eq.(^) leads to = T^* [^l], which is analogous to the Orr-Sommerfeld equation 

in the incompressible shear flow. (See Ref. |^^.) The term ( ^5| ) effectively gives the 'viscosity' in the viscous fluid. 
The viscosity is important for the stability of the flows. For the small Reynolds number, the laminar flow is realized, 
whereas the flow becomes turbulence above a critical Reynolds number. The effective Reynolds number in Eq.(|45|) 
is estimated from dimensional argument as 

Re-^. (46) 

This is roughly the square of the ratio of the g-mode frequency to rotational one. The viscosity term will play a key 
role on the singular point of the first-order equation, but the consequence is not clear at moment. It is necessary to 
explore further how the effective Reynolds number should operate in the stability and so on. 



VI. CONCLUDING REMARKS 



In this paper, we have explored an effective theory to describe the axial oscillation on a slowly rotating stars. The 
approximate equation governing the oscillation is constructed from the Einstein equations. The equation is derived by 
assuming that the angular dependence of the oscillation is dominated by a single component of spherical harmonics. 
The assumption in general breaks down as the evolution of oscillation. There are coupling terms of order e in the 
rotating fluids. Other oscillation patterns with different spherical harmonics will gradually be produced through the 
rotational coupling. For this reason, the equation is valid for small t, and can be used to examine the time evolution 
as the initial-value problem. The rotational effects up to third order are involved in this paper, so that the regime of 
application is enlarged. The equation derived here is also irrelevant to the singular point found in the first order one. 

The equation also shows a remarkable property. It is evident that the axial oscillation strongly couples to g-mode 
oscillations. Viscosity-like term arises from the polar pieces related to the g-mode oscillations. The 'viscosity' term 
originated from considering the sub-system only, i.e., a single component of the spherical harmonics. There are no 
production and extinction in the whole system, but, e.g., 'energy' of a component is partially transfered to the others. 
This transportation is regarded as dissipative effect so far as a particular sub-system is concerned. The term also has 
a significant implication. The condition A — Q in & good approximation for cold neutron stars, so that the coupling 
may be neglected. On the other hand, it is not clear that the condition holds, in particular for newly born hot neutron 
stars, in which the r-mode instability sets in. 

In this paper, we concentrated the equations only inside a star, to be more precise, the equations for the region 
A ^ 0. The pulsation equation derived here should be solved with appropriate boundary conditions. The boundary 
conditions are determined from matching with the equations outside, or regularity conditions. For example, the 
regularity condition of a function $ is given by $ r'~^ near the center. Depending on details of the stellar structure, 
the solution for the A ^ region may be matched to the solution for isentropic region, A = 0. Furthermore, the 
interior solution should be matched to the exterior one at the stellar surface. The exterior perturbation equation in 
vacuum is not derived here, but the form should be reduced to the wave equation describing gravitational wave. The 
perturbation equation should be solved by out-going boundary condition at infinity. One question may arise. Is it 
possible to calculate the radiation reaction at this order? Newtonian estimate indicates that the back-reaction is of 
order e2m+2 ^^j. ^ > 2. Our expansion of the rotational parameter is limited to the third-order, and higher order 
corrections are necessary to examine the effect in a consistent way. The radiation-reaction effect is also a kind of 
dissipative one, so that the accurate evolution for a long period is necessary. 
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APPENDIX A: EQUILIBRIUM CONFIGURATION OF A SLOWLY ROTATING PERFECT FLUID 



We here summarize the equiUbrium of a slowly rotating star to explain our notation. The equilibrium state 
with uniform angular velocity Q ^ 0{e) can be described by stationary and axisymmetric metric g^j,^, 4- velocity 
= (u*, 0, 0, u"^), pressure p and energy density p of the fluid. The rotational corrections up to 0{e^) is needed to 
assure the consistency in our analysis. The metric is given by 

2e' 



ds^ = -e^[l + 2{hQ + h2P2)]dr + e^[l + (mo + m2P2)\d 



W{1 + 2k2P2) < de'^ + sin^ 9 



dc^-\u + W^- W3^^ 1 dt 
smO dti 



(Al) 



where Pi = Pi{cos9) (/ = 2,3) is the Legendre polynomial of order I. The metric functions introduced above obey 
the following ordering in e: a; ~ 0(e), ho, h2,mo,m2, k2 ~ 0{e'^), Wi,W3 ~ 0{e^). These are functions of radial 
coordinate r only. The components of the 4-velocity are 



The pressure and energy density are respectively given by 

P = Po + {P20 + P22P2}, Pa + {P20 + P22P2}, 



(A2) 



(A3) 



where po and po are the pressure and energy density of non-rotating fluid. The centrifugal force of order alters the 
configuration shape, which corresponds to the quantities in the braces. The functions P20,P22, P20 and P22 are related 
with the metric functions of order e^. We rather use the metric functions to eliminate the pressure and density of 
order in the oscillation equations. 
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APPENDIX B: THE SECOND-ORDER TERMS 
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-167r(po + po)-^{4e'^wra7 - 2rv'{w + fi)^ - r^i/'ro'(ra7 + 3Q)} 



-u;^ - 2(r + -)e-''uw' - i^rr^e"^^'^ - 2/i + /2 + /. 



+ 



2 , 2 _ 19 2 . ,2 

24 



Z(i + 1) 



Z(/ + l) 

-A + /2 - 3/3 - A 



71 



24 



2e' 



2A 



+ Z(/+1) 3-TO0+ 1 



.2\ ,.,11 e^^o 



e a; — 



/in 



-IdTvp'oe^ ( e^mo + ye-'nj^^ + 2r^e-'vj'^ + /i - ^A + A, 



where 



(B2) 



A = 47r(p[, + 3p[,)(3e^6 + 2rh-''w^) - 16tt{pq + pofr^e^-^w^ - 12 (^1 

+^ (1 - 3e^ + 47r(po + Po)r^e^) 6 + - 3e^)e-''K72 + - r^e"^) e-'w' 



2 \ e 



in 



12 



+■ 



2v'e 



3 



/4 = I2 ( ^2 - -6 ) - 47r(po + Po) ( ( - + 4 ) tnV^e- - 



« ^12 

O 2A 



(B3) 



(B4) 



(B5) 
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3. Coefficient of J/f \ (3i 



Pi = C3 



- C2 



16 



A/2 /^A/2„2\' 



87re' 



2A 



(Po +Po) 



jy/j2j,4 



+ Ci 



"l2x 









- 18 



- 12 



Xfc2 



l{l + I)j2r2 



l{l+l) 



6^&6 + 12-^^^ 



I2OVK3 



j2r2 /(/ + I)j2r2 



X^'(c.r2)' 



+ 2x(e"^)' K + 



fw 1 /mo\'} 81^1+4814^3 327r 2 / 3 A 4 ,,3 



4. Coefflcient of /^z 

/32 = 4c2X [47r(/9o +Po) {4z<7(wr)' - vj'Lor} r'-^e^" - w' {uor)' j'^r'^] 



+ cixr^e " 



+ m 



-Y(Po + 3Po)t^V+n7'|(^^ 

X r 5 2 44 

47r(po + Po) S + y 

47r(p^, + 3p[,)xtJ7r4e-^ Qzi7 - ^x) 



5 \ /14 97 , ; X 
r 24)^+(y'^-24-^^'^ 



14 



6 



vjwr w vor w cor + „ , 

6 3 z/'^j^ 



3 19 2 

8^^" - 48^" 



9 1 

■2„2 J ^ /2„ . /„ 



XJ r <i -OT r - -X w - —X »^ 



r 11 8 1 1 

+47r(po +Po)xr-^e"'' Uzuryzur - —X'^r - -x^r + -j^xi'^e'")' - r^X^'r"^ + 



'„2„-A 



+ \ X-'^ + ^Xr e 



/in 



-(p^, + 3p^,)(4x-n7VVe-'' 



1 ,0,, 4 4 2 (w\' 7 i^' 2 / 4 J^V , 



+ 



2 1 
3-^ ~ 2 ' 



oX - 



wW2e"'' + -^{^''^'T - 8x'w + 4x'w'r)cc7j2r2. 



/33 = C2 



5. CoefHcient of /Js 



167rr2e-'^Xw(r-^w)'(p^, + 3p(,) - 167r(po +Po)e-''x 167r(po + Po)r-^e^n7(tu + Q) 



10 



+2r'^e^flzu ( 1 + 



e^po - 6ru){r^wy -^zu'iru'vcj + 20 - rw') 



+r^zu'{'lzu + brzu') I + — (r^e"'')'xwn7' - 4e"VxW(47r(po +po)r'^to + e"'^(ra;)') 



+ Cl 



(po +Po)0^ 1127rr^e^n7XPo - 4e^x'^ - 2r!^'(9n7 + w)x + 6rxro' + Tr^i^'x^' 



+r^i/'wX } + ^re'^ux [2w+ [1 



rv 



+ 



m 



167r, , 



-47r(po +Po)e t<7l — r + — r row H ^j"^^"^ 



^— (5ro + w) (l - ^ ) 
2 V rv' 



+ 



+r20ro'2e-" + 4xe-" (^e^ro^ - - 51 + 852 - 53 - (^1 - e^-'Ow^ 

l{l + - x)| (^1 - e^-''a;2 - ^moj + y (po +Po)r2e^-'^ro{n72 - A{w - oj)x} 



+ ^e"^{47r(87r(po +Po)r'x + X - ^)p[)r-' - x} ( + — j'ro^ 
+ i{(4 + u'r - 167rp^,r3)x + x'^} [ i^^o ] +K\-92+93 



J.2Q 

«(Z + 1)51 = 47r(p^, + 3p^,)<j ^(8ro + 30)6^^2 + ^e^-%w + n){Aw + 30) 



+^e-''r^vj\A{r^w)' - llOr) ^ + ^(po + Po)'e^-''r4ro2(8ro + 110) 
+47r(po +Po)e^|-|:(5ro - 3a;)6 + ^^'j'(8ro - 110)(12ro2 - (rro')^) 



Oe- 



,,/2 



- (-226^070) + 2ri/'(20ro2 + Uvuoj + 3^^) + 8rww' + (32tu + 13w)r2zy'ro') 



45 6 



16 



16 



+ 1 7 r'^ + (3e^ - 1 ; -n^^ i-^ f \ — re-'^Owro' 

rv' rv' 4:v' \ v' 



rv' 8 



+^^(23ro - 9uj)ojw' + ^^^(437k7 - 403a;)ro'^ + -^r^j^ro'^fllO - Sfrro)'), 
6 96 12 



11 



2r2 V 2^ / 6 12 



+4n{p'o + 3p'o)e^l^(2x + §) 6 + yin(r^w(0 - 4x) - 2x(r^a;)') 
_lr2g-.j^2(^ + f2)U - ^{po +pofr^e^-''w{Q.w - 4x{w + 20)) 



--e^x6 - ^^^e"" (12(07 + 2a;)tux + 4r{zu + 3uj)xzu' - {5zu + 2u>)r'^zu'x' 



r 3 

p-^ / 1Q 4 \ / , 47 n 

+^ - + 8x + ^(f^ - 4x) fe2 - ir^ e--- e-^n + + w 

\ 8 rv' J \ 8 ri^' 

+ ^e"''x's e"^z/'ro'(rro + (r^w)') - roro' + 2!y'tt7a; + ri^'w'{3w - w) L (Bll) 



-^jr- (Po + Po) ( 2(0 + ro)2 + rw'isn + zu) j 





4 2 




2^ 






- —UJ - 










v'r 







(B12) 



6. a[<i,if£'] 



where 



'(1) 



(1) 



'(1) 



1^^ 



Ai = Q± 
+S± 



2e' 



-u; + 647r(po+Po) (^■^j 



167re^(po +Po) <{ X - ^-72 (^^ + x) - ^ 



2e^ 
3 



1 :\oo 



Bi = Q±B3 + S± 



^^^ipo+Po)^2x{e^ +ru') + 2e^w-r'^u'm' - {2m + Lj) ^1 - ^ri^'^ rv' 



-(2ro + a;) 



-UJ 



2e^ ( rv' . 



+2xe^ (87rO(po+Po)-^) 



2xe^<^87rO(po+Po)- 1 



(B13) 



(B14) 



(B15) 



(B16) 
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Bo = 



S± 



-T± 



B3 = (2n7 + Lu) 



XUJ + 47r(po +Po) 



xBz - 2xe^ 



(72 ) +4^(Po+Po)^^ 



47re^(po +Po) 



— (l-e^) + — I w 



H (Po +P0) n7 



(B17) 



(B18) 
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